Universal dielectric loss in glass from simultaneous bias and microwave fields 
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We derive the ac dielectric loss in glasses due to resonant processes created by two-level systems 
and a swept electric field bias. It is shown that at sufficiently large ac fields and bias sweep rates the 
nonequilibrium loss tangent created by the two fields approaches a universal maximum determined 
by the bare linear dielectric permittivity. In addition this nonequilibrium loss tangent is derived for 
a range of bias sweep rates and ac amplitudes and show that the loss tangent creates a predicted 
loss function that can be understood in a Landau-Zener theory and which can be used to extract 
the TLS density, dipole moment, and relaxation rate. 
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Low temperature two level systems (TLS) in amor- 
phous solids have attracted growing attention recently 
due to their performance limiting effects in superconduct- 
ing qubits and resonators for quantum computing [3-0| 
and kinetic inductance photon detectors for astronomy 
[j]. TLS, commonly represented by atoms or groups of 
atoms tunneling between two configuration states (see 
Fig. [TJ 5]), are found to decrease the coherence of qubits 
due to microwave absorption. Amorphous dielectrics 
with TLS are commonly found in Josephson junction bar- 
riers and wiring crossovers @, Q , but TLS are even found 
on the surfaces of the resonators with coplanar supercon- 
ducting electrodes on crystalline substrates [4(. 

Despite theoretical and experimental studies of amor- 
phous solids over recent decades the identification 
of the tunneling systems and an understanding of their 
non-equilibrium phenomena have limited our ability to 
predict their influence on devices. In previous nonequi- 
librium studies the application of a dc bias voltage step 
with a small ac field frequency uj <C , results in a time 
dependent permittivity which indicates complex behav- 
ior of the TLS Initially after the electric field step, 
the permittivity increases quickly followed by a slow, log- 
arithmic, return to the equilibrium value. The observed 
behavior is interpreted as a consequence of the bias field 
induced change in the TLS density of states caused by 
long-range TLS interactions [1, [l(| ■ A numerical treat- 
ment of this low frequency phenomena compares the re- 
sult to interacting and non-interacting TLS theory 11 j. 

A recent experiment measures the high frequency 
(Huj 3> fcsT) ac loss of dielectric films with a simul- 
taneous parallel electric field bias. The study adjusts 
the sweep rate of the bias and the ac field amplitude 
such that they are able to probe various nonequilibiurm 
regimes of the films, and this provides an opportunity to 
understand nonequilibrium glass properties in the high 
frequency regime. The bias sweeps are created by apply- 
ing different dc voltage steps to the film through low-pass 
filters with time constants which are large compared to 
the resonator response time and the relaxation time of 



the TLSs. 

In high frequency steady state measurements of amor- 
phous dielectrics the resonant loss caused by TLS lowers 
(saturates) at large ac fields because the TLS Rabi fre- 
quency is larger than the TLS relaxation rate [H, [H [HI . 
However, in the nonequilibrium study with a large ac 
field the addition of a sufficiently large bias field sweep 
rate causes the loss tangent to increase back to the value 
found for the small ac field steady state [12 1 . 

In this paper we will show that the experimentally- 
observed maximum in nonequilibrium loss tangent should 
universally reach the value found for small ac fields in the 
steady state since it can be derived from the conventional 
tunneling model of amorphous solids. In addition we will 
derive the loss for a range of ac fields and bias rates and 
show that it follows a theory which depends on Landau- 
Zener transitions. The derivation of the nonequilibrium 
loss allows one to predict the TLS density, dipole moment 
and relaxation rate for a given experimental result. 




FIG. 1: The potential to a tunneling two-level system in an 
amorphous solid. A is the energy difference between the left 
and right well states when isolated, which are coupled with 
tunneling amplitude Ao. 

Each TLS can be characterized by its asymmetry, A 
and tunneling amplitude Ao [f|, which determine its en- 
ergy, E = \/A 2 + Ag. The TLS distribution is given 
by the universal law, P(A, A ) = j^, reflecting the ex- 
ponential sensitivity of the tunneling amplitude to the 
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two-well configuration. TLS interaction with an exter- 
nal electric field, F, is determined by its dipole moment 
p and contributes to the asymmetry energy A(F) = 
A(0) - 2pF. Following we assume the TLS dipole 
moment to be a randomly oriented vector with a fixed 
size p ~ 1 — 10D, which is assumed in the majority of 
theoretical models (e. g. [l3| )■ The use of different mod- 
els, e.g. a Gaussian distribution of dipoles will not change 
the main result, i.e. the universal value of the loss tan- 
gent obtained at large ac fields and bias rates. The loss 
tangent obtained here as a function of the bias rate and 
ac field amplitude depends on the dipole size, such that 
the loss tangent obtained in this work can be generalized 
for a distribution of dipoles. The interaction of TLS with 
the environment, i.e. phonons and other TLS, are char- 
acterized by its relaxation time T± and coherence time T 2 

In the absence of the external dc bias field the inter- 
action of TLS with the ac field F ac results in a field en- 
ergy absorption which can be described by the loss tan- 
gent defined as tan(<5) = p-, where e = e' + ie" is the 
complex permittivity. We consider the regime of a very 
high ac field frequency, to >• ^f- > j-,^-, respectively, 
which takes place for a typical experimental frequency 
v ~ 5GHz (jp ~ 200mK) and a very low temperature 
T < 30mK (2S3| ? considered in this letter. In this regime 
relaxation and decoherence for resonant TLS are deter- 
mined by spontaneous decays of excited states, so that 
T 2 « 2Ti Q. 

The spontaneous emission via photons can modify the 
ordinary decay via phonons in a given experiment, and 
both will contribute to relaxation. However, the cou- 
pling of the thermal bath is assumed to be sufficiently 
high such that the TLS far from resonance are initially 
in the ground state. As described below, this results in a 
TLS that will irreversibly (spontaneously) emit if coher- 
ently excited after crossing through resonance with the 
ac field, such that the coherent deexcitation processes at 
the crossing can be ignored. 

The contribution of a given TLS to the dielectric losses 
is determined by its interaction with the external mi- 
crowave field, F ac cos(wt), and environment, which are 
defined by the TLS Rabi frequency Qr ~ and re- 
laxation time T\. In the resonant regime, with Or <C u> 
and negligible ac Stark effects [lj|, one can express the 
TLS interaction Hamilonian h using the rotating frame 
approximation as 

h = -nn R s x — [e — hu)s z , n R = fw os (60^, 
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rain and maximum 
= E and the latter 



Here the minimum relaxation time Ti 
Rabi frequency £Iro are defined at Ao 
also uses a zero for the angle between the dipole p and 
the ac field F nr . 



Solving the Bloch equations for the TLS density matrix 
11] and averaging the result over the TLS distribution 
P(A, A ) = results in the following expression for the 
loss tangent [E II O 



tan(J) 



4ir 2 P oP 2 tanh(^v) 
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where where e is the relative permittivity (permittivity 
in cgs units, in our case of interest T 2 ^ m i n = 2Ti im j n ). 
Here the effective Rabi frequency after averaging over the 
angle 9 is Qm ^£l fl0 - The high field saturation of the 
loss tangent is caused by a TLS Rabi oscillations which 
are faster than the relaxation to the environment. The 
temperature dependence of the steady state loss reflects 
the thermal dependence of the TLS population difference. 

For a swept bias field parallel to the ac field the 
energy of the TLS depends on time as E(t) — 
a/ (A - 2pF bl as(t)) 2 + Af. A TLS contributes to the mi- 
crowave absorption near the resonance E w fvjj. Then 
one can approximately represent the energy of a TLS in 
the resonant form 



E(t) = huj + hv(t - t ), 



A 



v=vo Y-\j^) cos(0) ' hvo = 2p — (3) 

where to defines the time when the exact resonance takes 
place. Two resonances are generally possible for each 
TLS A - 2pF bias (t) = ±^{huj) 2 - Ag. We assume that 
the time between two resonant passages ^ is longer than 
the TLS relaxation time 7i. m i„ <wso that we can treat 
them independently. 

We begin the consideration of the non-equilibrium loss 
tangent with the oversimplified case when the relaxation 
and decoherence rates are small during resonance passage 
and can be neglected. Then the loss will take place after 
the resonance crossing events induced by the ac and bias 
fields (see Fig. [2]). In this regime TLS can be described by 
the wave function amplitudes in the ground and excited 
states (ci,C2), respectively. The modified wave function 
(ai,a 2 ) = (cie 1 "'/ 2 , c 2 e~ lult / 2 ) taken within the rotating 
frame approximation Eq. ([1} satisfies equations, 



da\ ,v(t — to) .Qr 
dt 2 2 21 

dd2 ,v(t — to) .flu 
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dt 2 2' 
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which are equivalent to Landau-Zener transition dynam- 
ics of a two-level quantum system (see Fig. [2 [l7|). If 
at t = — oo only the ground state is populated | a\ \ 2 = 1, 
then after the level crossing, t = oo, one has 
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The 'imaginary part of TLS dielectric constant e is de- 
fined by the reactive component of the TLS dipole mo- 
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FIG. 2: TLS energy spectrum as a function of time, induced 
by the bias field application. The ground and excited states 
are coupled by one photon transitions, described within the 
rotating wave approximation. 



merit, pi OS s, which can be expressed as 

a*(t) ai (t)). 



Plos 



-i^p(al(t)a 2 (t) 



(5) 



This expression should be averaged over TLS parame- 
ters including their energies, E — h{uj + v(t — to)), tun- 
neling amplitudes Ao and dipoles p. The integration 
over energy can be performed analytically employing the 
fact that (oJ(t)o 2 (i) - a* 2 (t) ai (t)) = if^Q ax | 2 - | 
a 2 | 2 ). Then this integration of Eq. <j5j) is equiva- 
lent to the integration over time, dE = vdt. It results 

in the Landau-Zener change in population difference, (g]). Indeed, at the crossover v 
2H " 



The intermediate regime, vq ~ ^fj , can be studied only 
numerically. The result of a numerical calculation of the 
Landau Zener formula of Eq. ([6]) is shown in Fig. [3] for 
dielectric losses as a function of the dimensionless sweep 
rate £ = 2 J?>° . This result is shown with the correspond- 
ing steady state non-linear loss for r^oIi.miraT^rmn = 81 
by a black line at 0.127. The two calculations do not 
agree in the asymptotic slow sweep limit, because the 
Landau-Zener calculations ignore relaxation and deco- 
herence processes during resonant passage. 

Relaxation and decoherence processes must affect di- 
electric losses at small field sweep rates [l!|. In fact in 
the case of fast relaxation f2j™Ti im i n T 2 , m i n <C 1 one can 
expect that the linear regime result will be valid at all 
bias field sweep rates (here the time T2,min stands for the 
inverse maximum decoherence rate including all possible 
decoherence channels though the results are presented 
for the low temperature limit T 2 w 2Ti). In the oppo- 
site, strongly nonlinear limit the microwave absorption 
should be saturated for small sweep rates, vq, and should 
collapse to the equilibrium case of Eq. @. In the steady 
state equilibrium (vq = 0) the non-linear microwave ab- 
sorption comes from the energy domain | E — huj |< 

hQuQ^J^^- If during the time T l min the change 

of TLS energy due to bias field sweep, SE ~ hvoTi. m i n , is 
small compared to the size of the domain, then one can 
ignore the field sweep and use the equilibrium result, Eq. 
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The average TLS loss tangent 
can then be written in terms of the previous expression 
integrated over the distribution in tunneling amplitudes 



tan((5) — 
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, Uro the equilib- 

rium non-linear loss tangent given by Eq. ^ and the 
non-equilibrium loss tangent given by Eq. (|7|) become 
equal to each other within the accuracy of a logarith- 
mic factor, on the order of unity. Thus one can quali- 
tatively approximate the dielectric loss behavior for dif- 
ferent Landau-Zener parameters, expressed through the 
dimensionless field sweep rate £ — 

l 



where < ... > indicates an average taken over the dipole 
direction with respect to the field. 

In the limit of a large bias sweep rate, J7|jq «o,one 
can approximate the exponent in Eq. (j6]) as e~ x ~ 1 — 
x. Then all integrals can be performed analytically such 
that we obtain the identical result as the steady state 
linear response limit (F ac = 0) of Eq. ([2]), tan(<5) = 
47r . This universal result is a consequence of the 
Fermi Golden Rule in a linear response limit, which does 
not depend on the nature of the 8- function broadening, 
determined by either the decoherence rate or the 

energy sweep rate (h^/vo). 

In opposite (adiabatic) regime, vq << one can es_ 
timate the integral in Eq. (j6]) with logarithmic accuracy, 
meaning that dx^—^, — ~ hi(a), as 



lty parameter rj 



and nonlinear- 
by either using the 



tan(M = ^^^ln( e - 1/4 ^ 
e irn 2 R0 \ 2u 
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steady state limit, Eq. ([2]), for small sweep rates, £ -^i ij, 
or Landau-Zener relaxation-free limit, Eq. ([S]). The lat- 
ter regime can be characterized by the asymptotic be- 
havior for £ <C 1 (Eqs. (JT))) or by the large bias rate limit 
in the opposite case (see Fig. [3]). 

We calculated the non-equilibrium dielectric losses us- 
ing a full numerical solution to the Bloch equations for 
each TLS [SEE!]. Monte-Carlo integration of the re- 
sults over the TLS distribution are shown by the thick 
blue line in Fig. [3] The result is consistent with the 
predicted behavior in the fast sweep rate limit and the 
steady state results of the slow sweep rate limit. 

Using this theory one can experimentally extract the 
dipole moment p, the TLS density Pq, and their relax- 
ation time, Ti, separately, in the low temperature limit 
under consideration (T 2 ~ 2Ti). Experiments can cre- 
ate a known bias sweep, ji" 3 , and ac field F ac and find 
the TLS dipole moment p that correctly sets the Lan- 
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dau Zener parameter, £, to agree with Fig. |3l The di- 

p 2 

mensionless parameter, — — , can be found independently 
from the loss tangent measurements either in the intrin- 
sic equilibrium regime, or in the strong nonequilibrium 
limit found for £ >> 1, which with the above informa- 
tion of p allows one to separately find Pq and T\ . Recent 
experiments with microwave ac and field bias contain the 
needed experimental regimes £ = 10~ 6 — 10 2 and r\ 3> 1. 
A full comparison of these experiments to this theory will 
be completed in a separate work. 

The theory is restricted to the low temperature limit 
where fcgT <C fku. At higher temperatures and low bias 
rate Uo^i.nun ^ the result remains applicable after 
multiplying the losses by the thermal population differ- 
ence factor tanh ( 2 fcsr ) ■ ^ larger bias rates the popu- 
lation difference cannot equilibrate to the instantaneous 
energy E(t), and the loss is determined by an earlier (and 
higher) energy E(t—T{). As a result the non-equilibrium 
loss tangent can exceed its steady state linear response 
limit if the TLSs are already thermally excited. 




FIG. 3: The non-equilibrium loss tangent as a function of 
the dimensionless field sweep rate, £ = , in the case 

^Ko2i,min22,uiin = 81. The thick blue line shows the Monte- 
Carlo averaged solution of density matrix equations. The 
solid black line shows the asymptotic behavior which includes 
the Landau-Zener solution (Eq. (JSJl), with the constant loss 
resulting from steady state behavior. The dashed lines ex- 
tend Landau-Zener and steady state behaviors to the whole 
domain of sweep rates. 



It is interesting that in the adiabatic regime many TLS 
are brought into their excited state creating a remarkable 
population inversion as in the rapid adiabatic passage 
regime leading to the phonon enhancement [2G . |21| . 

In conclusion we propose a theory to explain the effect 
of a sweeping electric field bias on the ac resonant loss in 
an amorphous dielectric. If the field sweep rate is very 
fast the loss tangent reaches a universal value even in the 
strongly non-linear regime of high amplitude ac fields, 
in agreement with recent experimental observations 12 1 . 
At slower bias sweep rates a strongly non-linear regime 
takes place, in which the loss tangent decreases with de- 



creasing the sweep rate until the saturation at the steady 
state non-linear limit Eq. ([2]). This nonequilibrium loss 
tangent theory can be used to interpret and extract many 
TLS properties in various contexts and provide a better 
physical insight into them. 
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